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ABSTRACT 

We have identified the presence of large-scale, low-frequency dynamo cycles in a long-duration, 
global, magnetohydrodynamic (MHD) simulation of black hole accretion. Such cycles had been seen 
previously in local shearing box simulations, but we discuss their evolution over 1,500 inner disk 
orbits of a global 7r/4 disk wedge spanning two orders of magnitude in radius and seven scale heights 
in elevation above/below the disk midplane. The observed cycles manifest themselves as oscillations 
in azimuthal magnetic field occupying a region that extends into a low-density corona several scale 
heights above the disk. The cycle frequencies are ten to twenty times lower than the local orbital 
frequency, making them potentially interesting sources of low-frequency variability when scaled to 
real astrophysical systems. Furthermore, power spectra derived from the full time series reveal that 
the cycles manifest themselves at discrete, narrow-band frequencies that often share power across 
broad radial ranges. We explore possible connections between these simulated cycles and observed 
low-frequency quasi-periodic oscillations (LFQPOs) in galactic black hole binary systems, finding that 
dynamo cycles have the appropriate frequencies and are located in a spatial region associated with 
X-ray emission in real systems. Derived observational proxies, however, fail to feature peaks with 
RMS amplitudes comparable to LFQPO observations, suggesting that further theoretical work and 
more sophisticated simulations will be required to form a complete theory of dynamo-driven LFQPOs. 
Nonetheless, this work clearly illustrates that global MHD dynamos exhibit quasi-periodic behavior 
on timescales much longer than those derived from test particle considerations. 

Subject headings: accretion, accretion disks — black hole physics — magnetohydrodynamics (MHD) 
— X-rays: binaries 



1. INTRODUCTION 

The standard physical model of black hole disk accre- 
tion accounts for the transport of angular momentum 
through correlations in magnetohydrodynamic (MHD) 
disk turbule nce driven by the magne torotational insta- 
bility (MRI, |Balbus fc Hawley|[T99l] ). While the linear 
behavior of this instability is analytically tractable and 
some analytic progress has been made in evaluating its 
saturation behaviors (e.g., iLatter et al. 1120091 iVishniac 



2009| |Pessah fc Goodm an 2009; Pcssah 2010), numen- 
cal studies are crucial tor understanding the nonlinear 
evolution of the MRI. Local (i.e., shearing box) simula- 
tions of weakly mag netized accretion, first conducted by 



Hawley et al. 



(1995), have been instrumental in facilitat- 
ing numerical studies of the MRI by permitting smaller 
domains and larger timesteps at a given resolution than 
their global counterparts. Global disk simulations, on 
the other hand, have been invaluable in elucidating the 
macroscopic aspects of accretion since they incorporate 
more natural boundary conditions and allow large-scale 
conservation behaviors and development of radial struc- 
ture. A crucial question is in what regimes do local sim- 
ulations serve as a true microcosm for global disk be- 
haviors. If, for example, large-scale magnetic structures 
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are important in disk coron ae, a s has been suggested by 



Blackman fc Pessah| (|2009j) and |Beckwith et al. (|2009j, 
one might worry that real accretion d i sks are intrinsically 

(2010) have shown 



non-local. Likewise, Sorathia et al. 
letic li 



that global magnetic linkages are important even though 
the evolution of fluid stresses in subdomains of a global 
simulation are well represented by local simulations with 
net magnetic flux. So while the character of MRI turbu- 
lence in global simulations appears to be correctly cap- 
tured by local simulations that include stratification and 
net magnetic flux, it remains unclear how behaviors in- 
volving the development of large-scale field correlations 
will translate from one simulation regime to the other. 

We report in this paper on a long-duration global sim- 
ulation of black hole accretion that confirms the exis- 
tence of an interesting phenomenon that previously had 
only been seen in local disk simulations. Specifically, we 
have detected prominent low-frequency "dynamo cycles" 
in the azimuthal magnetic field evolution of a simulated 
global accretion disk. Dynamo cycles are commonly in- 
voked as the explanation f or the observed 22-year so- 

Parker|[l955| |Babcock||1961 
w by Baliunas & Vaughan 



lar magnetic cy cle (see e.g. 
Lcighto^| JT969 



or the review by [ 



1985[ ). More directly relevant to our work, however, is 
the fact that such cycles have appea red in many shear- 
ing box simulations of accretion disks ( Brandenbur g et alT| 
r995l|Stoneet al.|1996[|Miller fc Stone|2TOH'lurner|iJ004[ 
Hirose et al.||2006[ |johansen et al.||2009| JSuzuki fc Inut-| 
suka|2009HShi et al'|2010||Gressel|2010[|Davis et al.|2010| 
Simon et al.||2010j ). These simulated cycles are typically 



seen to have periods on the order of tens of local or- 
bital periods, with the exact number varying somewhat 



2 



with the details of the simulation (the vertical domain, in 
particular, was a limiting factor in many of the ea rliest 
simulations). As |Shi et all (|2010|, iGressell (120101), and 



Davis et al. (2010) discuss, the generic "butterfly pat 



tern can be attributed to the vertical rising of azimuthal 
field due to the combined effects of dynamo action near 
the disk midplane and the Parker instability at higher 
elevation. Until now, however, such features have not 
been reported in global simulations. 

In analogy to the results of local simulations, our global 
disk simulation produces dynamo cycles with oscillation 
frequencies that are ten to twenty times lower than the 
local orbital frequencies for a large range of radii. Fur- 
thermore, the radial extent of our simulation captures the 
sharing of dynamo power at peak frequencies across rel- 
atively large radial ranges. As such, dynamo cycles pro- 
vide a tantalizing source of variability at frequencies com- 
parable to astronomically observed low-frequency quasi- 
periodic oscillations (LFQPOs) in multiple galactic black 
hole binary candidates. 

We proceed with a discussion of our numerical model 
(§2), followed by a detailed description of the resulting 
dynamo cycles (§3). We then compare our global simu- 
lation to published local simulations and discuss broader 
observational implications (§4), followed by our conclu- 
sions (§5). 

2. NUMERICAL MODEL 

Our fully three-dimensional magnetohydrodynamic 
(MHD) simulation employs a modified version of the 
p ublicly available ZEUS-MP (version 2) code, described 



in |Stone fe Norman| (fl992a|) , |Stone fc Norman] ( |1992 b[), 
andlHayes et al.|(|2006D. ZETJS-MP uses an Eulerian 



nite difference scheme accurate to second order in space 
to solve the equations of ideal compressible MHD, 
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We employ a gamma-law (7 = 5/3) gas equation of state. 
Timesteps are set by the usual Courant condition, and 
a protection routine prevents the density and pressure 
from reaching artificially small and/or negative values. 
The only adjustments we have made to the fundamental 
ZEUS-MP algorithm involve this protection routine and, 
as described below, the introduction of modified gravity 
and the gas cooling function, A. 

Gravity is modified in our simulation to emulate the 
relevant effects of gen eral relativity through a pseudo- 
Newtonian potential ( Paczynski fc Wiita||1980[ ) of the 
form: 

GM GM 
* = -p3W (6) 



This approach accurately captures for a Schwarzschild 
spacetime the position of the innermost stable circular 
orbit (ISCO) at r = 6r g , the period of which is tisco ~ 
61.6 GM/c 3 in this potential. 

We initialize our computational grid using spherical 
coordinates (P,6>, 0) that span P g [4r g , 400r g ], 9 £ 
[0.057T, 0.957t], 4> € [0, 7r/4). The grid is non-uniform, log- 
arithmically increasing in R with a maximum R reso- 
lution of AP = 0.025 r g at the inner edge of the grid. 
The zone aspect ratio is approximately AP : RA9 : 
RA(f> = 3:1:6 everywhere within seven scale heights 
above/below the disk midplane (i.e., at 6 = n/2), and 
each scale height is resolved with 25 computational zones 
in this region. Outside of this region, the 9 resolution 
logarithmically increases outward. The total grid size is 
N R x N 9 x = 512 x 384 x 64 = 1.26 x 10 7 zones. 
Standard ZEUS-MP boundary treatments are employed, 
with a restricted boundary condition that permits out- 
flow only in the ±P directions. Reflecting conditions are 
used near the coordinate pole in 9 (as in De Villiers & 
Hawley 2003 for example), and periodic conditions are 
applied in <p. 

The initial conditions correspond to a thin, axisym- 
metric disk of constant midplane density and radially 
decreasing pressure: 



p(R, 9) = po exp 



cos 2 9 



2(h/r) 2 sin- 2 



and 



p(R,( 



GMR(h/r) 2 sin 2 9 



(P 



2r g ) 2 



p(R, 



(7) 



(8) 



where po is the initial density in the disk midplane and 
h is the effective scale height. The disk aspect ratio 
is initialized to h/r — 0.05 everywhere, and a cooling 
function A is implemented to maintain this aspect ra- 
tio with a cooling time (r coo i) that is related to the lo- 
cal orbital period (r or b) using an ap proach similar to 
that described in Noble et al. (2009). The exact form 
of the cooling function is A = /(e — etargVicooi, where 
/ = 0.5[(e — etarg)/|e — etargl + 1] is a threshold function 
that enables cooling only when the internal energy e is 
greater than the target energy etarg (and which is equal 
to zero when etarg > e). The target energy is chosen so 
that e talg oc pw 2 (/i/r) 2 , which comes from the assump- 
tion that c s ~ (h/r)vrf, in thin disks. Estimating the cool- 
ing time as the thermal timescale of the disk, we choose 
Tcooi = Torb/a = 1 0r or b, which corresponds to a Shakura 
fe Sunyaevf (|1973| alpha disk with a = 0.1. While the 



true effective alpha parameter varies spatially and in time 
over the evolution of an MHD disk, this choice provides 
an adequate order-of-magnitude estimate for the imple- 
mentation of cooling. Additionally, we ran a short test 
simulation that revealed that the frequency range of the 
dynamo cycle signal discussed in the following section 
was insensitive to the presence or absence of cooling, al- 
though the signal itself was more pronounced in the case 
with cooling. 

The initial velocity profile is entirely azimuthal and is 
set such that the effective centrifugal force balances the 
gravity of the central object in the disk midplane. The 
initial magnetic field is completely poloidal in orienta- 
tion and consists of a series of magnetic field loops that 



3 




2x10 



4X10 6x10 
Time (GM/c 3 ) 



8x10 



2x10 



4x10 6x10 
Time (GM/c 3 ) 



3x10' 



2x10 4x10 6x10 8x10 
Time (GM/c 3 ) 

(right). The field has been 



Fig. 1. — Space-time evolution of azimuthal magnetic field at R = 15 r g (left), 20 r g (center), and 25 r g 
averaged only in the azimuthal direction. Bright (yellow) regions indicate strong, positive fields while dark (violet) regions are strongly 
negative. The disk midplane is located at 8 = n/2 and A8 = 0.05 is one scaleheight. Time is shown both in units of GM/c 3 and in terms 
of the local orbital period, r or b. These "butterfly diagrams" illustrate that the azimuthal field changes sign over timescales on the order of 
ten times the local orbital period. 



span seve ral local scale heights in both height and width 
(e.g., see Reynolds k. Miller||2009[ ) . The average ratio of 
gas-to-magnetic pressure is initialized to f3 1000. 



3. RESULTS 



The following analysis treats the disk only after it 
has relaxed away from its initial conditions. Specifi- 
cally, we allow the disk to evolve for 150 ISCO orbits 
(> 9200 GM/c 3 ) and define t = to correspond to the 
end of this initialization period. The simulation is then 
followed post-initialization for more than 1500 ISCO or- 
bits, from t = to t k 9.8 x 10 4 GM/c 3 . During and after 
the initialization, the disk naturally evolves to a turbu- 
lent state as a result of the MRI. Deferring a detailed 
discussion of the full disk evolution to a future work, we 
focus here on an interesting set of coherent behaviors 
that emerge from this turbulence. 

3.1. Azimuthal Magnetic Field Oscillations 

Figure 1 shows the azimuthal field strength (also az- 
imuthally averaged) as a function of both time and el- 
evation from the midplane for three distinct radii. At 
all radii, the disk midplane is located at 8 — ir/2, while 
A6 = 0.05 corresponds to one scaleheight for fixed h/r. 
Each panel of the figure shows that the azimuthal field 
at a given location reverses sign multiple times. This 
variation can be seen for several scale heights above and 
below the disk midplane, forming a "butterfly" pattern 
analogous to that which has been observed in shearing 
box simulations of accretion disks. The period of field 
reversal is generally longer for larger distances from the 
central object, although there are also some common fea- 
tures seen at all radii. While it is difficult to tell from 
these diagrams whether the variability is truly periodic, 
it appears that the reversals in field orientation take place 
on timescales on the order of tens of local orbital periods. 
Near the maxima of the field cycles, the azimuthal field 
strength can be twice as strong as the RMS value of the 
azimuthal field measured over the simulation duration in 
the same region. 

To explore this variability in more detail, we show in 
Figure 2 the power spectral density (PSD), defined as 
P(v) — t]\f(y)\ , where 77 is a normalization constant 



and f(v) is the Fourier transform 

/» = / f{t)e- 2 ™ l dt, 



(9) 



of a time series fit). In each panel, the azimuthal mag- 
netic field has been azimuthally averaged and summed 
over a radial range comparable to the local scale height 
before the PSD is computed for a location two scale 
heights above the disk midplane. In this case, the time 
series is taken to be the duration of the simulation after 
initialization. All three panels show strong power en- 
hancements at frequencies roughly comparable to ten to 
twenty local orbital periods (i.e., 0.05 — 0.1 v Q xb, where 
t'orb is the local orbital frequency). It is interesting to 
note that all radii also show multiple peaks, suggesting 
that the phenomenon is not always simply related to the 
local orbital period. In fact, there is some overlap be- 
tween adjacent radii, as seen in the shared frequency of 
the strongest peaks for both R = 15 r g and R = 20 r g . 

As a coarse estimate of the significance of these peaks, 
we can compare their strengths to the mean power as 
estimated from nearby frequencies. For a PSD of a sin- 
gle time series, the mean is comparable to the standard 
devia tion of the power distribution 



(e.g. 



Press et al. 



1992|, implying that the ratio of peak-to-mean power 
can be used to estimate peak significance. In the case of 
R = 15 v„, for example, the two highest peaks are ap- 
proximately 15-20 times the mean power as extrapolated 
from nearby frequency ranges. The case of R = 20 v g 
is even more convincing as the primary peak is approx- 
imately 90 times stronger than the mean while the sec- 
ondary peak is roughly 40 times above the mean. Even 
if our approximations somehow underestimate the mean 
power by factors of a few, each radius features multiple 
peaks that stand significantly above the noise. 

To further evaluate the significance of these features, 
we also construct the average power spectral density 

N 

(PSD), defined as P(v) = (l/N) ^P;(z/) over a set of 

i=l 

N independent time series fi(t). This approach has the 
advantage of reducing the standard deviation of the PSD 
features at the expense o f the available frequency domain 



and resolution (see e.g., |van der Klisj[1989[ |Press et al. 
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Fig. 2. — Frequency-weighted PSDs (vP, in arbitrary units) of azimuthal magnetic field strength taken over the duration of the simulation 
at at R = 15 r g (left), 20 r g (center), and 25 r g (right), measured at two scale heights above the disk midplane and averaged both azimuthally 
and over a radial domain comparable to a scale height. Each radius features strong, multi-peaked signals corresponding to the oscillations 
seen in Figure 1. 
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Fig. 3. — Frequency-weighted PSDs (vP, in arbitrary units) of azimuthal magnetic field strength as in Figure 2, but constructed from 
averages over four independent time series. Each radius features strong signals suggesting that the oscillations seen in Figures 1 and 2 are 
robust. 



1992 Vaughan et al.|2003 |. Given the frequencies of in- 
terest aS3 - tne T toIaT^aiiration of the simulation, we can 
afford to take only N — 4 independent time series, which 
increases the significance of the averaged PSD by a fac- 
tor of two over the unaveraged case. Figure 3 shows the 
PSDs of the azimuthal magnetic field over the same re- 
gions depicted in Figure 2. As in Figure 2, we see peaks 
(now broadened, but at higher significance) that sit ap- 
proximately ten to twenty times below the local orbital 
period. While the peaks are sufficiently broad that they 
overlap for different radii, the frequency resolution of the 
PSD is insufficient to determine how significant this over- 
lap is. 

Figure 4 shows the PSD of the azimuthal magnetic field 
as a function of both frequency and radius, so that we 
are better able to evaluate the radial dependence of the 
power profile. As in Figure 2, the PSD has been com- 
puted over the entire time series. The vertical features 
in Figure 4 illustrate that power is shared at distinct fre- 
quencies across large radial intervals, with single peaks 
often stretching across radial ranges of 10 r g or more. 
Taken in aggregate, however, the power distribution re- 
flects the radial run of the orbital frequency. Specifically, 
the power is bounded by ~ ^orb/6 on the high-frequency 
end and ~ i/ or b/30 on the low. So even though a given 
power peak may radially span multiple orbital frequen- 
cies, the range of peak frequencies remains approximately 
proportional to the local orbital frequency. This pattern 
only stands out clearly from the noise for R 10 r g . In- 
ward of this region, the broadband noise associated with 
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Fig. 4. — PSD (P, in arbitrary units) of azimuthal field strength 
over a range of radii, as measured at two scale heights above the 
disk midplane and averaged azimuthally. Also shown as a dashed 
line is the local orbital frequency. That the dark (red) band of 
enhanced power runs parallel to the dashed line demonstrates that 
the range of oscillation frequencies is roughly proportional to the 
local orbital frequency. A given power peak, however, can often be 
seen to stretch across a large radial range. 



accretion across the ISCO masks any obvious trend. 

To further elucidate the nature of this variability, Fig- 
ure 5 presents analyses of the azimuthal field at R = 20 r g 
after various cuts and segregations have been imposed. 
In the left panel of Figure 5, we show the PSDs for this 
region when it is divided into two independent time se- 
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panel illustrates the time variation of the signal by comparing PSDs taken over two independent halves of the total simulation time (after 
startup). The middle panel shows PSDs corresponding to two distinct regions (two scale heights) above and below the disk midplane. The 
right panel shows PSDs for three and four scale heights above the midplane, which are similar in frequency to one another and to the signal 
at two scale heights. 



ries (Tl represents the first half of the simulation, T2 the 
second), each of which is approximately 5 x 10 4 GM/c 3 
in duration. Note that the peaks are not constant in 
frequency, suggesting that the multiple peaks in Figures 
2-4 are at least partly caused by frequencies migrating in 
time. It is tempting to claim that the peaks move from 
higher to lower frequencies over time, but it is challeng- 
ing in practice to follow individual peaks without a much 
longer time baseline. 

The middle panel of Figure 5 shows the PSD for the full 
time series, now comparing the regions above and below 
the midplane. These peaks, too, fail to perfectly align 
even though the disk starts from an approximately sym- 
metric state. This is perhaps not surprising given that 
features in the disk turbulence are also seen to evolve 
asymmetrically about the midplane. That said, this top- 
bottom asymmetry suggests that the exact frequency of 
a single field oscillation may be less useful as a diagnostic 
tool than the range of frequencies observed. The right- 
most panel of Figure 5 compares the PSDs for regions 
three and four scale heights above the disk midplane. 
The similarity of these peaks to one another and to the 
locations (if not the relative amplitudes) of the analogous 
peaks in Figure 2 suggests that, as expected, the variabil- 
ity on a given side of the disk features similar frequencies 
at different heights from the disk midplane. 

3.2. Observational Proxies 

Thus far, we have focused only on the behavior of 
the azimuthal magnetic field. We would also like to 
explore whether any derived quantities and/or observa- 
tional proxies oscillate in a similar manner. Figure 6 
shows PSDs of three derived quantities: the integrated 
R — 4> stresses 

W R ^ = J ( P v R v^ - B K B^/47r)dV, (10) 

the integrated Ohmic dissipation 

p„^/^2!* (id 

(where a is assumed constant), and the mass accretion 
rate 

M = - [ P v R dA R . (12) 

-'ISCO 



Integration ranges are provided in the caption to Fig- 
ure 6, and the RMS-normalized PSD for each quantity 
is computed over the entire simulated time series. Of 
these three quantities, only the total stresses feature os- 
cillations that stand above the local noise at frequencies 
comparable to those at which the azimuthal field oscil- 
lates. This is not completely surprising since the domi- 
nant Maxwell stresses depend upon the azimuthal field, 
but it is nonetheless interesting that a quantity spatially 
integrated over a large radial range, five scale heights, 
and the entire azimuthal range of the grid still selects a 
set of distinct frequencies. In contrast, both the Ohmic 
dissipation and accretion rate are dominated by noise, 
particularly at low frequencies. In practice, we would 
have to convincingly model and subtract what appears 
to be a red (i.e., Brownian) noise spectrum in the full 
data set to prove that any features in the full PSDs were 
significant. 

As was done in Figure 3, we also construct the PSDs 
for these three observational proxies as averaged over four 
independent time series to evaluate the significance of the 
peaks in Figure 6. For each proxy, we see features near 
the frequencies of interest (i.e., a few xlO~ 4 c 3 /GM), 
but they are in all cases comparable (within approxi- 
mately a factor of two) to the noise levels in the nearby 
lowest frequency bins, as seen in Figure 7. This is in stark 
contrast to the case of the azimuthal magnetic field (Fig- 
ures 2 and 3) , where a clear signal stood above the noise 
in both the full PSDs and the PSDs constructed from 
the subdivided series. Without better frequency resolu- 
tion from a longer simulation, it is difficult to conclude 
that any peaks in the observational proxies are present 
at high levels of significance. 

4. DISCUSSION 

4.1. Interpretation and Comparison with Local 
Simulations of Accretion Disks 

We propose that the magnetic field oscillations clearly 
seen in Figures 1-5 are the manifestation of a dynamo cy- 
cle phenomenon. As discussed in the Introduction, such 
cycles have been seen in local simulations of magnetized 
accretion that reflect a large variety of numerical algo- 
rithms and disk parameters. Our work has shown that 
similar cycles are indeed present in global simulations 
of black hole accretion disks. As in the local case, the 
azimuthal magnetic field cycles with frequencies approx- 
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Fig. 6. — Frequency-weighted PSDs (vP, where P is RMS-normalized) of the total stresses (left), Ohmic dissipation (center), and accretion 
rate across the ISCO (right). The stress and dissipation are integrated over a range that spans R = 15 — 25 r g , stretching in elevation 
from the disk midplane to five scale heights above it. The accretion rate is measured over a range that spans five scale heights above the 
midplane. While the integrated stress shows signs of the azimuthal field signals seen in Figures 2-3, the dissipation and accretion rate are 
noise-dominated . 
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Fig. 7. — Frequency-weighted PSDs (vP, where P is RMS-normalized) of the total stresses (left), Ohmic dissipation (center), and accretion 
rate across the ISCO (right), as in Figure 6 but averaged over four independent time series. While there are features in each panel near 
the frequencies of interest (i.e., a few x!0~ 4 c 3 /GM), these features are not particularly well separated from the noise. 



imatcly ten to twenty times lower than the local orbital 
frequency. Individual peaks in these cycles can share 
power across large radial ranges, while the range of cy- 
cle frequencies remains roughly proportional to the local 
orbital frequency. Additionally, these features persist for 
the full duration of the simulation (nearly 10 5 GM/c 3 ). 
This is much longer than the time expec ted to erase all 



memo ry of the initial field conditions (see Sorathia et al 
2010 for example) and is comparable to the radial drift 



timescale at r = 10 r g , suggesting that this behavior can 
be maintained even over a period during which the disk 
evolves significantly. 

Linking the amplitudes and frequencies of the observed 
dynamo cycles to the parameters of our global disk sim- 
ulation is, in practice, qui t e cha llenging. For example, 



Blackman & Brandenburg ( 2002 ) illustrate with a series 
ot numerical experiments that the resulting dynamo cycle 
frequency depends on the time evolution of the effective 
magnetic diffusivity. In a grid-based simulation such as 
ours, this will generally vary with both the grid resolution 
and flow details. Additionally, our simulated disk does 
not maintain a perfect steady-state. For example, the 
surface density at a given radius can vary by 50% during 
the course of the simulation as accretion proceeds. This 
variation is not necessarily monotonic, however, and is 
accompanied by variation in both the accretion rate and 
effective "alpha" parameter. It is therefore difficult to as- 
sociate the potential frequency drift seen in the first panel 
of Figure 5, for example, with any obvious trend in the 
disk, although a longer time baseline would potentially 
enable such an identification. As such, we take a phe- 
nomenological approach to the observed dynamo cycles 
and, having noted their similarity to a well-established 



aspect of local simulations, now discuss how they relate 
to astrophysical observables. 

4.2. Relevance to Astrophysical Black Holes 

The astronomical phenomena most obviously analo- 
gous to our simulated dynamo cycles are the LFQPOs 
detected in multiple black hole binary systems by X- 
ray satellites such as Ginga and the Rossi X-ray Tim- 
ing Explore r. As summarized in |Remillard fc: McClin-| 
tock (2006), LFQPOs appear in the nonthermal X-ray 



spectrum and have frequencies that range from 0.1 to 30 
Hz. Although the peaks are often seen to migrate quite 
rapidly in frequency, a given peak typically has a qual- 
ity factor (Q = v/Av) > 10. Observed LFQPOs can 
also be very strong, exhibiting RMS amplitudes of tens 
of pcrccnts. 

Our simulated dynamo cycles have three properties in 
common with observed LFQPOs. First, they occur in 
the expected frequency range. To choose an example, an 
observed 4 Hz LFQPO for a 10 M black hole corre- 
sponds to a frequency in natural units of v ~ 2 x 10~ 4 
GM/c 3 , which is representative of the range that we see 
in the simulations. Second, several of the simulated dy- 
namo cycle peaks have high quality factors. In practice, 
the duration of our simulation permits a maximum qual- 
ity factor (3 max ~ 5 at a frequency of 2 x 10~ 4 GM/c 3 , 
but this value is achieved at multiple radii in Figure 2, 
for example. Naturally, the quality factors of the peaks 
in the time-averaged PSDs shown in Figure 3 are lower 
(specifically, of order unity) as a result of the reduced 
frequency resolution available in this mode of analysis 
and, potentially, the non-stationarity of the time series. 
Third, the dynamo cycles are seen to occupy the magne- 
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tized, low-density region that has been traditionall y iden- 
tified with the hard X-ray emitting corona (e.g., Miller 
& Stone 2000). Thus, there is no difficulty in linking 
dynamo cycles to the region where we think observed 
LFQPOs originate. That noted, we should point out 
that our simulation does not include any treatment of 
realistic radiative processes. Neither do any of our con- 
structed observational proxies reflect the magnetic field 
cycles with a high degree of significance, nor do the prox- 
ies feature RMS amplitudes comparable to those seen in 
observed systems. As such, we cannot make any detailed 
predictions concerning how these dynamo cycles would 
manifest themselves observationally. 

Whatever the advantages and disadvantages of dy- 
namo cycles as a model for LFQPOs, it is worth briefly 
distinguishing them from the most popular alternative 
models of LFQPO production. First, that dynamo cy- 
cles are most prominent in the azimuthal magnetic field 
off the disk midplane shows that they ha ve none of the 
obvious characterist ics of trapped waves ||Katol|l9 90t or 
diskoseismic modes flNowak fc Wagoner 1991 , 1992, 1993; 

[2009| lO'Neill et ai.||2009| FH5otF^f 
est themselves m the hydrodynamic 



Reynolds fe Miller 
which should mani 
variables. Lense-Thirring pre cession has been explored 
in t he con text of LFQPOs by Ipser ( 1996 ) and Stella et 
al. (1999), for example , in the test particle limit and on 
a more global scale by Ingram et al. (2009), but our nu- 
merical model of a non-rotating black hole is obviously 
incapable of producing such effects. (In fact, given the 
observed dynamo cycles' locations and nature, one ex- 
pects a similar outcome from a fully general relativistic 
simulatio n.) The truncated disk mo del for LFQPOs pro- 
posed by Giannios & Spruit (2004) is also quite distinct 
from ours since dynamo cycles naturally produce oscil- 
lation frequencies much lower than the local orbital fre- 
quency even at tens of gravitational radii, requiring no 
disk cutoff. Finally, we note that dynamo cycles are not 
simply a subtle manifestation of the A ccretion- Ejection 
Instability (AEI, |Tagger fc Pellat|1999|) since the AEI oc- 
curs near the inner edge of the disk while our cycles are 
seen to originate from all radii that have had sufficiently 
many orbital periods over which to evolve. 

5. CONCLUSIONS 

We have described a global, numerical, MHD study 
of black hole accretion that has revealed an interesting 
pattern of magnetic variability that, until now, had only 
been seen in local shearing box simulations. The most 
important results of our study are summarized here: 

1) We have identified for the first time the presence 
of dynamo cycles in global simulations of black hole ac- 



cretion disks. These cycles manifest themselves as oscil- 
lations in the azimuthal magnetic field in a region that 
stretches from a few to several scaleheights above the 
disk midplane in elevation. Individual peaks in these cy- 
cles share power radially, while their frequency range at 
a given radius is found to be approximately ten to twenty 
times lower than the local orbital frequency. 

2) While dynamo cycles are easily seen in the azimuthal 
magnetic field, detecting cyclic variation in derived quan- 
tities is much more challenging. The integrated stresses 
feature variability that we have identified with the az- 
imuthal field behavior, reflecting the fact that distinct 
peaks in the PSD are manifested by a wide range of radii. 
The mass accretion rate and integrated Ohmic dissipa- 
tion, however, remain noise-dominated, and none of the 
three proxies we have examined feature RMS amplitudes 
as large as those seen observationally. 

3) We have discussed potential links between dynamo 
cycles and observed LFQPOs in black hole binaries. 
While dynamo cycles naturally produce oscillations at 
the appropriate frequencies and locations expected for 
LFQPOs, any complete theory of LFQPOs will need to 
address how the field oscillations generate an observable 
signature, how their periodicity varies in time, and, ul- 
timately, how low-frequency variability is connected to 
the inferred black hole accretion state. Additionally, the 
duration of our simulation limits our LFQPO peaks to 
have quality factors Q < 5 when estimated from a single 
time series. Such quality factors are still a factor of two 
below what is observed, implying that direct compar- 
isons of this nature with observed LFQPOs will require 
extended simulations that run for several thousand or- 
bits (or longer, if one wishes to average multiple PSDs 
together). Nonetheless, our work shows that dynamo cy- 
cles in global simulations can produce oscillations with 
frequencies much lower than any natural frequency in the 
test particle regime. 
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